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Abstract: A new signal processing algorithm for absolute temperature measurement using 
white light interferometry has been proposed and investigated theoretically. The proposed 
algorithm determines the phase delay of an interferometer with very high precision (« one 
fringe) by identifying the zero order fringe peak of cross-correlation of two fringe scans of 
white light interferometer. The algorithm features cross-correlation of interferometer fringe 
scans, hypothesis testing and fine tuning. The hypothesis test looks for a zero order fringe 
peak candidate about which the cross-correlation is symmetric minimizing the uncertainty 
of mis-identification. Fine tuning provides the proposed algorithm with high precision sub- 
sample resolution phase delay estimation capability. The shot noise limited performance of 
the proposed algorithm has been analyzed using computer simulations. Root-mean-square 
(RMS) phase error of the estimated zero order fringe peak has been calculated for the 
changes of three different parameters (SNR, fringe scan sample rate, coherence length of 
light source). Computer simulations showed that the proposed signal processing algorithm 
identified the zero order fringe peak with a miss rate of 3 x 10" 4 at 31 dB SNR and the 
extrapolated miss rate at 35 dB was 3 x 10" 8 . Also, at 35 dB SNR, RMS phase error less 
than 10" 3 fringe was obtained. The proposed signal processing algorithm uses a software 
approach that is potentially inexpensive, simple and fast. 

Keywords: Phase delay estimation, white light interferometry, signal processing algorithm, 
fine tuning 
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1. Introduction 



Although fiber optic interferometric sensors offer the possibility of performing measurements with 
very high sensitivity and resolution [1], they suffer from problems such as complex signal processing 
techniques, quadrature point stabilization, and uncertainty as to whether an increase or decrease in the 
value of measurands has occurred. In order to fully utilize the capability of fiber optic sensors, a 
different sensing principle termed as "white light interferometer" or "white light interferometry (WLI)" 
was developed [2]. 



Figure 1. All Fiber White Light Interferometer. 
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From the beginning of this research, an all fiber white light interferometry (AFWLI) absolute 
temperature measurement system as shown in Figure 1 was selected as the application area of the 
proposed signal processing algorithm. White light interferometry departs from the conventional 
interferometry in that it uses a broadband light source. SLD in Figure 1 represents a superluminescent 
diode (SLD) used as a broadband light source and PD1, PD2 denote photodetectors 1 and 2, 
respectively. In the AFWLI two fiber Fabry-Perot interferometers (FFPI, sensing FFPI and reference 
FFPI in Figure 1) and a Mach-Zehnder Interferometer as a processing interferometer (scanning 
interferometer in Figure 1, hereafter termed as MZI), which has two piezoelectric transducers 
(hereafter termed as PZT) in its two arms, are connected in tandem. The sensing and the reference 
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interferometer output signals from PD1 (Equation 1 for PD1) and PD2 (Equation 2 for PD2) are given 
respectively by: 



7 s (<D P ) = l + |exp 



7 fl (O p ) = l + |exp 
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where Lc is the coherence length of light source, Op is the OPD (Optical Path-length Difference) 
between two arms of Mach-Zehnder processing interferometer and Os, Op are the round trip phase 
shifts for the sensing and reference FFPIs, respectively. In Equation 1 and Equation 2 it is assumed that 
light source has a Gaussian power spectrum. A constant d.c. voltage Vdc (100-150 volt, Figure 1) is 
applied to the PZT1 in one arm to coarsely match the OPD of the MZI to that of the sensing FFPI. And 
an alternating ramp voltage Vsaw (Figure 1) is applied to PZT2 in the other arm to scan the processing 
interferometer so that OPD Op between two arms of MZI can be varied over a certain range. 



Figure 2. Output of sensing and reference FFPI from AFWLI. 
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This AFWLI for temperature measurement produces two fringe scans, one from the sensing FFPI 
and another one from the reference FFPI, as shown Figure 2. In AFWLI, the sensing FFPI is exposed 
to the temperature Is to be measured and the reference FFPI is protected from environmental 
disturbances but exposed to the known reference temperature T R . Now, assume that the known 
temperature of the sensing FFPI and the reference FFPI are Ts and T R , respectively. When the phase Op 
of MZI is scanned and exactly matched to that of the sensing FFPI (the reference FFPI), a zero order 
fringe peak of the sensing FFPI (the reference FFPI) is produced at certain Op = Ops (Op = Op«), as 
shown in Figure 2. In Figure 2, Ops (Op«) denotes the phase of the processing interferometer 
producing the zero order fringe peak of the sensing FFPI (the reference FFPI). If we can identify the 
phase difference O rf = Ops - Opp (this is possible by the proposed signal processing algorithm which 
will be shown later in this article) then O rf is mapped to the temperature Is and absolute temperature 



Sensors 2008, 5 



7612 



measurement is possible. This problem has been known as "Time Delay Estimation (TDE) [3] (or 
Phase Delay Estimation)". In this article, a new signal processing algorithm to estimate the phase delay 
O rf of AFWLI is proposed. This article consists of five sections. Section 2 describes the previous 
related works for time delay estimation methods. In Section 3, a signal processing algorithm to 
measure the phase difference <b d of AFWLI is proposed. In Section 4 the performance of the proposed 
signal processing algorithm has been demonstrated using computer simulations. Section 5 shows the 
comparison to the previous literature and Section 6 is the conclusion of this article. 



2. Previous works 



Two major classes of signal processing algorithms for WLI are the hardware approach and the 
software approach. Both approaches have a more or less "tracking zero order fringe peak" feature. 
Gerges proposed a hardware approach which locks the zero order fringe position of interferometer by a 
feedback loop [4]. An improvement of the sensitivity up to 1/240 fringe was claimed. To the author's 
best knowledge this method, while still dependent on the incremental characteristic of laser 
interferometry and not fully taking advantage of WLI's potential to identify the interference fringe, 
demonstrated the feasibility of locking and tracking the fringe peak for absolute measurement for the 
first time. 

There are many software algorithms to estimate phase delay O rf [3, 5]. Among many detection 
methods, the cross-correlation method dominates the field of phase delay estimation in practice due to 
its easier implementation [6]. Many other TDE methods are based on this algorithm. The cross- 
correlation method cross-correlates the two fringe signals is{n) and z«(n), which are sampled versions 
of Isi^p) and Ir{Q>p) respectively, into i (n) and considers the sample number argument n=rid that 
corresponds to the maximum peak in cross-correlation i (n) as the estimated time delay [7]. 

While WLI has the potential to identify the interference fringe order from the output pattern of an 
interferometer [8], it is difficult to distinguish the zero order fringe peak from its adjacent first order 
fringe peaks when noise is present in the interferometer output, as shown in Figure 3. From Equations 
1 and 2, the amplitude difference AI between the zero order fringe peak and the first order fringe peaks 
is represented as: 

f ^ Y 



AI = 



1 



1 -exp 



(3a) 



Clearly, if a system has a noise level which is equal or greater than AI the zero order fringe peak 
cannot be identified directly, simply through inspection of its amplitude. If the normalized zero order 
fringe peak value is defined as unit signal, a minimum signal-to-noise ratio SNR mm required to identify 
the zero order fringe peak [9] is given by 

" ( 9 V 

= -201og< 
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Representative values of the coherence length of different light sources like white light lamp, light 
emitting diode (LED), superluminescent diode (SLD), are about 10, 20, and 40 in terms of optical 
fringes. The SNR mm required to identify the zero order fringe peak by amplitude difference AI is given 
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from Equation 3 as 28 dB, 40 dB, and 52 dB respectively [9]. One disadvantage related to cross- 
correlation is the broadening of envelope (or peak) (hereafter termed as envelope) from Lq to -Ji L c 

due to the cross-correlation of two almost identical Gaussian envelope signals [17]. This results in 
higher SNR mm in Equation 3 requiring 6 dB more than before cross correlation. Then, with this 
broadening, the attainable resolution is often not better than one fringe and a precise peak location may 
be somewhat questionable. This difficulty has inhibited the application of fiber optic sensors using 
WLI, for example, absolute OPD determination [9]. 



Figure 3. Output of white light interferometer ( Is in Equation 1 or Ir in Equation 2). 
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Additionally, once the zero order fringe peak is identified, then for a more accurate sub-sample 
resolution time delay estimation we will have to use interpolation which is possible by either quadratic 
interpolation in time domain or frequency domain zero-padding [10, 11]. But quadratic interpolation 
uses three cross-correlation values centered at the estimated peak of cross-correlation. This method has 
a shortcoming of bias caused by time domain sample rate [12, 13] and the difference between true 
peak shape and the fitted quadratic polynomial. In frequency domain zero-padding, the number of 
zeros are padded in the middle of Fourier Transform of cross-correlation i(n). Notice that zero-padding 
in frequency domain increases the discrete frequency by a certain factor which eventually results in 
time-domain interpolation in cross-correlation signal i(n). Zero-padding in frequency domain is a 
useful tool to improve the peak location accuracy, but it increases the computational complexity [14] 
and storage requirement associated with inverse Fast Fourier Transform (FFT) operations [15]. 

Notwithstanding the above mentioned shortcomings, cross-correlation is a still useful tool for time 
delay estimation as shown that cross-correlation with no pre-filtering is an optimal maximum 
likelihood estimator to estimate the time delay between two similar signals if the noise processes wr 



Sensors 2008, 5 



7614 



(n), ws (n) of signal in (n), is (ri) are white noises and at least one of signals has high signal-to-noise 
ratio greater than 10 dB [16]. 

The outcome of the time delay estimation depends on the combined performance of coarse 
estimation (zero order fringe peak identification) and sub-sample resolution estimation of time delay. 
In this article, a new signal processing algorithm which can accurately identify the zero order fringe 
peak of cross-correlation i (n) of two fringe scan output signals of a WLI is proposed. This algorithm 
still uses a cross-correlation technique taking advantage of simple implementation. But this algorithm 
combines the hypothesis test as a coarse estimation to reduce the possibility of mis-identification of 
peak with fine tuning algorithm as a sub-sample resolution peak estimation to overcome the 
shortcomings of quadratic interpolation or frequency domain zero-padding. 

The proposed signal processing algorithm uses a software approach, which is potentially 
inexpensive, simple and fast. And, this proposed signal processing algorithm has a low peak mis- 
identification rate of 3 x 1 0" 4 at 3 1 dB SNR and has a high precision fine tuning capability down to 5 x 
10" 4 fringe as will be shown from the computer simulation results. 



3. Proposed Signal Processing Algorithm 



The proposed signal processing algorithm consists of five steps applied to sampled signal of WLI 
fringe scans. They are: 

1) Normalization and cross-correlation 

2) Peak and zero crossing detection 

3) Matched filtering 

4) Hypothesis test 

5) Fine Tuning 

Each procedure is explained briefly below. 



3.1. Normalization and cross-correlation 



As a preliminary procedure, the output of photodetector signals Is (®p) and Ir (0,p) are sampled and 
normalized respectively into fringe scan i s (n) and 

i s (n) = exp<j - : — - >cos^.7r(n-n s )/ f s }+w s (n) (4) 



7rL c I A 



and: 
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cos{2;r(n -n R )/f s }+ w R (n) 



(5) 



where ns is zero order fringe peak sample number for sensing FFPI, n« is zero order fringe peak 
sample number for reference FFPI, f s is the sample rate in samples per fringe (or samples/fringe). A 
normalization procedure is carried out by removing the d.c. component (constant "1" in Equation 1 
and Equation 2) of each fringe scan. ws(n) and Wr(h) are the white noise related to the fringe scan is{n) 
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and i R (n) respectively with zero mean and variance u w . The phase delay between the sensing FFPI 
and the reference FFPI is defined in terms of samples as: 

n d= n s~ n R (6) 
After normalization, is (n) and in (n) are cross-correlated into i(n) and normalized again. Cross- 
correlation i(n) can be expressed in mathematical form with its zero order fringe peak p 0 a.t n=rid as 
shown in Equation 7. 

00 

*(»)= Z**(*)*s(« + *) 

k= 

2n(n-n d )lf< 



exp 



cos{2x(n - n d ) / f s }+w(n) (7) 



In Equation 7 the effective coherence length Lc, e ff is given as yflL c [17] and w(n) is a white noise 
of i(n) . Exponential term in Equation 7 is termed as "the envelope" of cross-correlation i(n). This 

cross-correlation improves SNR at zero order fringe peak up to 14 dB [17]. One disadvantage related to 
cross-correlation is the broadening of peak from Lq to -JlL c due to the cross-correlation of two 

(almost) identical Gaussian envelope signals [17]. For example, this results in higher SNR m m in 
Equation 6 requiring 6 [dB] more than before cross-correlation when Lc=26X, but this will be 
compensated by 14 dB SNR improvement at the zero order fringe peak [17]. After cross-correlation, 
the task of the signal processing algorithm becomes to find a zero order fringe peak po=rid (global 
maximum) of cross-correlation i (n) correctly. 

3.2. Peak and Zero-crossing Detection 

At this stage, all the peaks p,, all the minima qi and all the zero crossings z, of the cross-correlation 
i(n) are detected and registered. Every peak is labeled as p t where its value po=nd is zero order fringe 
peak position, p.\ is negative first order fringe peak position, p\ is positive first order fringe peak 
position in terms of sample number and so on. Also every zero crossing between peak p { and minimum 
qi is detected and labeled as z,- where j= 1,2,3... in terms of sample number. Then zero crossing position 
Zj with sub-sample resolution is calculated [9] by using interpolation formula as shown in Equation 8, 
Equation 9 and Figure 4. 

Zj =r+ ' (r) for \i(r)\ > \i(r + 1)| (8) 

and 

J i(r + \)-i(r) 1 1 1 1 

Linear interpolation is the reasonable method because the cosine function crosses the zero 
essentially as a straight line as shown in Figure 4. 
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Figure 4. Linear Interpolation and zero crossing. 
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Table 1 shows one example of peak and zero crossing table. Note that in this table peak position is a 
integer number but zero crossing is not necessarily a integer number due to the interpolation 
calculation shown in Equation 8 and Equation 9. 



Table 1. Example of peak and zero crossing table (fs =20 [samples/fringe]). 
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Figure 5. Least Square Fit of zero crossing positions. 
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Then the zero crossing period b is calculated by fitting all the zj into the form of: 

y(n) = bn + a forn = 1,2,3,. ..N z (10a) 

by applying least square fitting to the distance between all pairs of zero (z/ and z j+1 ) as shown Figure 5 
where Nz is the total number of zero crossings in i(n). In Figure 5 the slope b of the y in) is the 
estimated number of samples inside half fringe (or half the sample rate) and the estimated sample rate 
is given by: 

f s = [2b] (10b) 

where [2b] is the closest integer number of 2b. 



3.3. Matched Filtering 



A matched filter is the optimum filter to maximize the SNR of signal out of the matched filter in the 
presence of additive stochastic noise at input signal [18]. Consider the case where the input signal is 
s(t) and n{t) is the white noise with zero mean, variance g 1 ^ related to s(t). In this case, matched filter 

theory states that the maximum SNR at the output will occur when the matched filter has an impulse 
response h(f)=s(to-f) that is equal to the time-reversed version of the signal waveform s(f) to which 
matched filter is matched and time-delayed by certain to [second]. When input signal s(f) is time- 
limited signal existing only when 0< t <T then the matched filter is defined by [18]: 

h(t) = s(T-t) (11a) 

Note that the output so(t) of the matched filter, which is the convolution of s(t) and h(f)=s(T-t), is 
given by: 

(•00 (• 00 

s o (t)=h(t)*s(t)=\ h(A)s(t-A)dA = \ s(T - X) s(t - X)dX (lib) 

JO JO 



Figure 6. Matched filtering. 
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In Equation 1 lb so(t) is the cross-correlation of the time-reversed version of input signal s(t) (i.e. s(- 
t)) with h{t)=s{T-t), which is denoted as R s (. t )h(t)(f). But s(-f) and h(t) are identical two waveforms with 
time distance T separated and so(t)=R s (. t )h(t)(t) becomes R ss (t-T) when R ss (t) denotes the autocorrelation 
R ss (t) of input signal s(f). Then so(t) has its maximum value when t=T and this maximum value is same 
as the R ss (0) of autocorrelation R ss (f) of input signal s(t). 

Thus if we process a signal-plus-noise with a matched filter, the largest peak outputs due to the 
signal will correspond to R ss (0) and the SNR [18] of the largest peak output of matched filter is given 
by: 
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SNR 



output of matched filter 



(12) 



If we define the SNR of input signal s(t) as the ratio between the largest value a of input signal s(t) 
and oat, then SNR improvement due to the matched filtering is given as: 



SNR 
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(13) 

'AT °V a 

At this stage, following the concept of matched filter, only one fringe of cross-correlation signal i 
(n) (i.e. signal between i(p,) and /(/?,+/) of i (n)) is considered as the signal to be matched and one 



period of cosine function, cos 



' 2m 
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is considered as the matched filter, 



i M (n) = cos 
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n = 0,l,-f s 



(14) 



Then, for a given one fringe of cross-correlation signal i{n), signal out of the matched filter is 
maximized when the time delay between two signals is f s and maximum value is calculated as in 

Equation 15. 

n=p, n=p i \ JS J 

for given J; 's 0=0, ±1, ±2,...). The envelope function of i{n) is slowly changing in the vicinity of the 
zero order fringe peak of i (n) and Equation 15 can be rewritten as: 
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(16a) 



where A, in Equation 16 is the average value of exponential envelope function between p t and p i+ i as 
shown in Equation 16b. 

i(Pi)+i(p i+ i) 



h = 



(16b) 



For the case of a zero order fringe peak i(p 0 ) of i(n), average value ho of envelope between p 0 and 
p\ is approximately "1" because z (p 0 ) = 1 and i (p 0 ) ~ i (pi). Then, Equation 13 becomes: 



S3VR 



imp, matched 




= 



(17) 



a 1 
and SNR of z'(/?o) is improved by the factor of: 

201og(V^)*5[dB] (18) 

In the above Equation 16a any J, value is the weighted integration of one fringe between two peaks 
of i{n) and also can be considered as maximum value out of matched filtering which has a improved 
SNR over i(p,) as shown Equation 13. So, J t values will be used in the hypothesis test instead of peak 
value i(p,) of i{n) as will be shown later (see Section 3.4.). One useful property of J, is Jo= J-i, J\= J-i 
. . . (Ji= J-i-i) due to the even symmetric property of i (n) about n=po= rid. 
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3.4. Hypothesis Test 

In hypothesis test, signal processing algorithm chooses the nine biggest peaks of i (n) as zero order 
fringe peak candidates. Ideally, p/s of/=0, ±1, ±2, ±3, ±4 are selected. Hypothesis test presumes that 
each candidate peak is the zero order fringe peak and calculates the parameter gipj): 

GO GO 

g{p j ) = Y J {Jj+i-J H -i)=Y J dU,i) for / = 0,±1,±2,±3,±4 (18) 

i=0 i=0 

where j of notation d(j,i) corresponds to the candidate peak pj on the hypothesis test and i is the 
distance from the zero order fringe peak candidate in terms of fringe. Then g(po) is expressed as: 

g(p 0 ) = d(0,0) + d(0,l) + d(0,2) + ... 
= (J 0 -J_ l )+(J 1 -J_ 2 ) + (J 2 -J_ 3 ) + ... (19) 

Ideally all the values of d(j=0,i) for the zero order fringe peak po is zero (Figure 7) due to the 
symmetric property of J,= J.i-i and the zero order fringe peak candidate producing | g(pj) | =0 is 
announced as the estimated zero order fringe peak p 0 . But, practically the zero order fringe peak 
candidate producing minimum | g(pj) | is announced as the estimated zero order fringe peak p 0 . Note 
that ideally zero order fringe peak po happens at n= rid, first positive order fringe peak p\ at n= rid+ fs 
and so on. 

Figure 7. Distribution d(j,i) of noise-free cross-correlation. 
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3.5. Fine Tuning 

Figure 8 shows the vicinity of zero order fringe peak of cross-correlation i(n). From Figure 8 it is 
clear that the discrete sample zero order fringe peak position pd of the cross-correlation i(n) identified 
by the hypothesis test is not necessarily same as the "true zero order fringe peak position", p t and the 
goal of fine tuning algorithm is to find the distance between pd and p t . 
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Figure 8. Illustration for true zero order fringe peak and discrete sample zero order fringe 
peak of i(n). 
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First, fine tuning algorithm assumes that the cross-correlation i (n) is represented as: 



i(n) = exp 



2n[n- p)l f s 

nL C,eff I A 



fcos- 



2 n(n-p) 

fs 



(20) 



Note that p t in Equation 20 is not necessarily a integer sample number. Also, fine tuning algorithm 
generates test cross-correlation, i test (n) defined by: 



W(«) = exp 



exp 



2n(n-p d )l f s 



2x(n-n d )lf s 



>cosj 



\2x(n-p d )\ 



>cosj 



fs 

\2x(n-n d ) 



fs 



(21) 



In Equation 21 it is assumed that the discrete sample zero order fringe peak position p d happens 
ideally at n= n d and also that the approximate value of the coherence length of the light source, L c is 

known priori to us. This is possible by counting the number of fringes within full width at which the 
intensity of interferogram using a particular light source decreases down to e x of its maximum value. 
And value of f s is given from f s = [lb] as shown in Section 3.2. 

Then, the distance between p d and p t can be found by calculating J(M) for varying integer number 

M: 
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J(M) = ^i(n)i test (n-MAn) 



(22) 



-co<M<co 



where i tes t (n-MAn) is a computer-generated test cross-correlation with its zero order fringe peak at p d 
+MAn and An is the desired fine tuning resolution in sample. 



An = 



N. 



[sample] 



(23) 



sub 



In Equation 23 N su b is the number of sub-divisions in one fringe (the desired fine tuning resolution 
in fringe). Then J(M) becomes: 

2n(n-p t )l f s 



J(M) = X exp 
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exp 
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2n{n-n d - MAn)] 

7s l 



(24) 



Interpretation of Equation 24 is that fine tuning algorithm generates its own test cross-correlation, 
itest {n) with the true zero order fringe peak positioned at nd +MAn and calculate J{M) for varying M. 
Then J(M) has a peak value at the certain Mf which makes M/An ~ p t - nd because the best similarity 
between i (n) and i test (n) is attained when nd +M/An ~ p t . Then, fine tuning algorithm announces {nd 
+M/An) as the estimated true zero order fringe peak location p t . 

Calculating J(M) for the range of-co<M<+co takes many calculations and is time-consuming. 
The true zero order fringe peak p t normally is within several samples of the discrete sample zero order 
fringe peak pd. Especially, when SNR of the fringe scans is high enough, the true zero order fringe 
peak p t is ideally within the half sample of the discrete zero order fringe peak pd. So, for the above case 



J(M) is calculated only for the range of 



N sub /f s 
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~N sub /f s ~ 


2 
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N sub /f s 



Equation 24 can be expressed as shown Equation 25 where | ^- Js 



for faster signal processing and 
is the closest integer number of 



AM) 
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= £ exp 
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jUn-ri) 
1 fs J 



exp 



2n{n — n d — MAn)/ f s 



-|2" 






>cos| 







2n{, 



n-n r 



-MAn)} 



fs 



(25) 



Additionally it must be emphasized that the principle of fine tuning in Equation (24) and Equation 
(25) can be used for matched filtering in order to obtain the best similarity between the signal inside 
the one fringe of cross-correlation (signal between any peak i(pi) and i(pi+\)) and the matched filter 
z'm (n) and to maximize the SNR improvement at peaks by matching the phase between the discrete 
sample z'-th order fringe peak z'(p ; ) and the very first value of the matched filer z'm(«)- 
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4. Computer Simulations 

The proposed signal processing algorithm was verified using computer simulations. To see the shot- 
noise limited performance of the proposed signal processing algorithm, the normalized AFWLI fringe 
scans, is(n) and i R (ri) were computer-generated using Equation 4 and Equation 5 respectively and shot 
noise was added to the AFWLI fringe scans, is(n) and i R (n) instead of white noise. In the computer 
simulations the zero order fringe peak position po,s (and po,s) for is{n) (and i R {n)) were randomly 
selected as real number. is{n) and i R {n) were cross-correlated into i{n). Then, p t is calculated as (po.s- 
Po,r) and zero order fringe peak pd= m of On) becomes the integer part of p t . The coherence length of 
is{n) and i R {n) were chosen as L c =261 to simulate the coherence length of the commercially available 
SLD like OKI OE350S from Oki semiconductor. For all the computer simulations presented in this 
article parameters are fixed as follows, unless the specified parameter is varied for a certain range and 
circumflex notation a on the top of the parameter denotes the estimated value of that particular 
parameter calculated by the proposed signal processing algorithm. 

• Sample rate fs=l6 [samples/fringe] 

• SNR=30 dB 

• Effective Coherence length Lc, e ff = V2x26 X « 36/1 

• Size of fine tuning step An=l/1000 [fringe] 

4. I. Simulation: Miss rate simulation 

In the first simulation miss rate (misidentification rate) of the proposed signal processing algorithm 
was tested at different shot noise levels. The SNR was varied from 1 dB to 28 dB with 1 dB separation 

and a set of 10000 simulations was carried out at different SNR. When the position difference between 
the estimated zero order fringe peak p d and the computer-generated zero order fringe peak p d is 

bigger than half fringe (8 samples), the zero order fringe peak is considered to have been misidentified. 
Table 2 shows the miss rate of the proposed signal processing algorithm and miss rate is defined as the 
ratio of number of miss to 10000. 

In Figure 9, miss rate was plotted along with Bit Error Rate (BER) of binary fiber optic 
communication system [19]. In Figure 9 (or Table 2) we have a rule of thumb that every dB 
improvement in SNR (over the range of 26 ~ 3 1 dB) produces approximately one order of magnitude 
improvement in error rate. This kind of behavior is also the case for binary fiber optic 
communication system (two orders of magnitude improvement in error rate for the binary fiber optic 
communication system). 

To extrapolate the miss rate beyond the range of 10" 4 on the BER curve, data points in abscissa in 
Figure 9 were left-shifted by -14.2 dB (by trial and error) until data points of miss rate between 16 dB 
~ 31 dB were visually fitted on the BER curve. Then four more data points were extrapolated beyond 
the data point of miss rate at 31 dB SNR on the BER curve as shown in Figure 10. It is predicted that 
miss rate will be 3 x 10" 8 at SNR of 35 dB. 
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Table 2. Miss rate of the proposed signal processing algorithm as a function of SNR. 
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Figure 9. Comparison between miss rate and BER. 
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Figure 10. Extrapolation of miss rate on BER curve. 
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4.2. Simulation: Root mean square (RMS) error of the zero order fringe peak identification 



After the zero order fringe peak was identified, fine tuning was calculated for resolution 
enhancement. Phase error 0 CTTOni between computer generated zero order fringe peak p t (or p 0 ) and 
fine-tuned zero order fringe peak p t (or p 0 ) was calculated. Phase error O errorii was averaged over 30 

simulations at a given SNR and this average gives out the root mean square (RMS) error of the zero 
order fringe peak identification. 
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Figure 11. Change of RMS error as a function of SNR (10-40 dB). 
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Figure 11 shows the change of RMS error along with the SNR in the range of 10 dB to 40 dB and 
Figure 12 shows the change of RMS error along with the SNR in the range of 30 dB to 40 dB. It is 
shown that the minimum SNR required to achieve RMS error less than 10" [fringe] (which is the fine 
tuning step size) must be greater than 35 dB SNR. 

Figure 12. Change of RMS error as a function of SNR (30-40 dB). 
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4.3. Simulation: High Precision Performance of Fine Tuning 



In this simulation, the SNR was fixed at 30 dB and RMS error was compared between two cases. 
One is the case where the estimated zero order fringe peak p 0 was fine-tuned for resolution 
enhancement and the other case was where the estimated zero order fringe peak p 0 was not fine-tuned. 

Figure 13. Comparison between RMS error with fine tuning and RMS error without fine 
tuning. 
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As shown in Figure 13, without fine tuning, RMS error of p 0 is approximately 03/fs fringe and 
RMS error of p 0 totally depends on the sample rate. But this RMS error was reduced down to -0.0015 
fringe when p 0 was fine tuned. Note that RMS error was not sensitive to the sample rate fs over the 
range between 10-20 samples/fringe. This feature predicts that signal processing algorithm can use 
lower sample rate fs for the faster signal processing, but still maintain lower RMS error by using fine 
tuning, which is comparable to the error of higher sample rate fs. In the proposed algorithm, the cross- 
correlation between is in) and in(n) is most time-consuming part of signal processing algorithm, but the 
sample rate fs can be lowered to speed up the signal processing without losing higher sample rate RMS 
error. One strong benefit of the proposed algorithm is that once zero order fringe is identified correctly 
then the resolution of the fine tuning can be adjusted as low as desired simply by changing the value of 
N S ub and high precision phase delay estimation can be obtained. 

4.4. Simulation: Effects of b ( f s ) 

In this simulation, estimated effective coherence length L Ce jj was set as 36 fringes, SNR=30 dB 

and estimation error of the zero crossing period b was varied from -10% to 10% when b=\0 (fs =20 
samples/fringe) to see the effect of estimation error of the zero crossing period b on the performance 
of the proposed signal processing algorithm. The estimation error of the zero crossing period b is 
defined by: 

Estimation error of b = - — — xl00[%] (27) 

b 

where b is the half the sample rate fs. In Figure 14, it can be shown that RMS error of fine tuning was 
not sensitive to the estimation error of b within the range of ± 6 % estimation error of b and 
especially, RMS error of 0.001 [fringe] was obtained within the range of 2% estimation error of b . But, 
on the range of±7%~±10% estimation error b RMS error increased dramatically to 0.25 fringe. 

Figure 14. Change of RME error as a function of estimation error in zero crossing period. 
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This is presumably due to the fact that i(n) and test cross-correlation i te st(n) (Equation 24 or 
Equation 25) started to be out of phase and RMS error increased fast. 

4.5. Simulation: Effects of L c 

This simulation is to show the effect of estimated coherence length L c on RMS error assuming that 
b = b , Lc=26X{Lc, e ff~ 36A) and SNR was fixed as 30 dB. Estimated coherence length L c was varied 
from 2\X to 30A, . In other words, estimated effective coherence length L c e <r was varied from 30>l to 

42A, and RMS error was observed. Figure 15 shows that RMS error turned out to be not sensitive to the 
estimation error in coherence length L c ej j . This is due to the property of cross-correlation. Cross- 
correlation is maximized when p t and pd +MAn are in phase as long as both cross-correlation i (n) and 
test cross-correlation i tes t(n) are symmetric. 

Figure 15. Change of RMS error as a function of estimation error in light source coherence 
length. 
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5. Comparison to the literature 



In this section a comparison to the previous literature data regarding the resolution of zero order 
fringe peak detection is given. Interestingly enough, reference [20] calculated the theoretical limit of 
scanning white light interferometry signal evaluation algorithm. In this reference the theoretical limit 

of resolution of the fringe order detection was given as: 

N 

Az„ * — (28) 



— f I (zfdz 
Az J - co 



1/2 



S t = S( Zi ) = l( Zi -z c ) + N{z t ),i=0, n (29) 

where Si=S(zi) consists of a fixed number of intensity values, typically taken at equidistant positions z ; , 
I(z) is the ideal input signal (i.e. fringe scan) to be subject to intensity noise N(z), z c of I(z) is the 
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position to be found using the evaluation algorithm (signal processing algorithm). And it is assumed 
that intensity noise N(z) is the constant noise value N over all samples. In our case, cross-correlation 
i{n) can be represented as the generalized equation as follows again: 



i(n) = 1 + 1 0 exp 
Then, the derivative of i{n) becomes: 
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( 2nn^ 
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and (z'(z)) is given as: 
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Then, substituting Equation 32 into Equation 28 produces: 
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(31) 
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(33) 



Io/N'm denominator in Equation 33 is the signal-to-noise ratio and Equation 33 becomes: 

jAzl ceff x 0.353x1.263 

Az c *^ ^ X (34) 

IJN 

As can be seen from the Equation 34, theoretical limit of fringe order detection is the function of 
sample rate, effective coherence length and SNR. Theoretical limit of resolution for the computer 
simulation shown in Figure 1 1 and Figure 12 can be calculated using the parameters of 16 sample per 
fringe (Az= X/\6), effective coherence length l ceff = 42 x 26 x A = 36 A . Then theoretical limit of 

resolution of zero order fringe detection in Figure 11 and Figure 12 is given as: 

V(A/16)(36/l) x 0.353 x 1.263 0.669x1 . t c 

Az c « m = (without fine tuning) (35) 

SNR SNR 

V(A/1000)(36A)x0.353xl.263 0.0846A , 

Az r « — — = (with fine tuning) (36) 

SNR SNR 



Figures 16 and 17 are the comparison of the computer simulation results and the theoretical limit of 
fringe order detection calculated by Equation 35 and Equation 36. As shown in Figure 16, computer 
simulation results produce a much bigger RMS error than the theoretical limit over low SNR range 
(10-30 [dB]). In comparison to the literature, it seems like that theoretical limit of reference [20] is 
more optimistic over low SNR range than the computer simulation performance. But, the proposed 
signal processing algorithm has proven to reach the theoretical limit of fringe order detection over the 
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higher SNR range (30-40 [dB]). This is not a surprising result because in higher SNR fine tuning 
algorithm will enhance the fringe order detection resolution down to the theoretical limit once zero 
order fringe peak is identified correctly. But, over low SNR range, fine tuning algorithm is not 
effective in enhancing the fringe order detection resolution because probably zero order fringe peak is 
misidentified and fine tuning is searched within the half sample of the misidentified zero order fringe 
peak. 



Figure 16. Comparison of computer simulation and theoretical limit (SNR: 10-40 [dB]). 
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Figure 17. Comparison of computer simulation and theoretical limit (SNR: 30-40 [dB]). 
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Increasing the fine tuning range will help to locate the zero order fringe peak correctly and lower 
the RMS error down to the theoretical limit, although this is time consuming. Optimum fine tuning 
range of M over low SNR range to find the correct zero order fringe peak will be the subject of further 
research. 
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6. Conclusions 

A new signal processing algorithm for white light interferometry has been proposed. The goal of 
signal processing algorithm was to find the time delay (phase shift) between two fringe scans which 
makes it possible to measure the absolute optical path length of a sensing interferometer. This new 
signal processing algorithm can be used for absolute temperature measure measurement by mapping 
the zero order fringe peak position of cross-correlation i(n) to the time delay between two fringe scan. 
Cross-correlation between two fringe scans utilizes all the photons in fringe scans effectively and the 
uncertainty of the zero order fringe peak mis-identification was reduced. Monte-Carlo simulations 
showed that the proposed signal processing algorithm identified the zero order fringe peak with a miss 

4 8 

rate of 3x10" at 31 dB shot noise and the extrapolated miss rate at 35 dB shot noise was 3x10" . Also 
resolution of less than 10" fringe was obtained at 35 dB shot noise (Figure 12). The fine tuning of 
signal processing algorithm requires some prior knowledge on the coherence length Lc of the SLD and 
sample rate fs. But the performance of proposed signal processing algorithm turned out to be not 
sensitive to the estimation error in L c and f s (within the range from -6% to 6% estimation error of 

b ). Also the proposed signal processing reached the theoretical limit of fringe order detection over the 
higher SNR range (30-40 [dB]). The proposed signal processing algorithm uses a software approach 
which is potentially inexpensive, simple and fast. As a whole, the proposed signal processing 
algorithm has proven to be a high precision signal processing algorithm for AFWLI phase (time) delay 
estimation. 
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Appendix A 



In this appendix, SNR improvement at the peak of cross-correlation i(t) is derived. First, z's(®) and 
ir (O) are given respectively as: 



where O p = 



ImLf 
X 



2mL 



< — 


<D P -O s 


2~ 

> 




tcL c 1 X 




< — 


~o P -o R ~ 


> 




7rL c IX 
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cos(O p -O s ) 
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(A-l) 
(A-2) 



X a ' K X 

Here L P , Ls, L R are the path length difference of the Mach-Zehnder processing interferometer, 
sensing FFPI (twice the Fabry-Perot cavity length), reference FFPI and n is the effective refractive 
index of core of the optical fiber (not to be confused with sample number n in i(n)). Notation 
<8>,* and <=> are reserved to denote the cross-correlation, convolution and Fourier transform 
relationship. 
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When L P is scanned L P is a function of time t and also Op is a function of time t. Then: 

<b p =2fft=w 0 t (A-3) 

where / (or wo) is a modulation frequency (or angular velocity) of fringe scan. Following the same 
token, (£>s, ®« are given by: 

® s =2xft s = w 0 t s (A-4) 
® R =2rft R = w 0 t R (A-5) 
and the relation between and tn is given by: 

t s =t R +t d (A-6) 

where ^ is the time delay between sensing FFPI fringe scan and reference FFPI fringe scan. 
Using the following substitutions: 

O P -O s =w 0 (t-t s ) = w 0 (t-t R -t d ) (A-7) 

O P -O R =w 0 (t-t R ) (A-8) 
L c = N f A (A-9) 

and considering the noise in is(®p) and /«($/>), time domain representation of the Equation A-l and 
Equation A-2 can be rewritten as: 

is(t) = e -^) 2 CO s(w 0 (t-t d )) + w s (t) (A-10) 

i R (t) = e- a ' 2 cos(w 0 (0) + w R (t) (A- 1 1 ) 

In Equation A-10 and A- 11 we used a substitution a = (wo/xNj) and it is assumed that 0^=0 
without loss of generality and that ws(t), WR(t) are white noise with zero mean and variance <j 2 n of is(t), 

iR(t) respectively and independent each other. Also, wo was set to "1" for normalization in order to 
make is(t) = e l at t = (N//2) x 2k. 

The first part of Appendix A shows that the increase of coherence length from L c to sj2L c due to 

cross-correlation. We know that cross-correlation and convolution are identical operations for the case 
where either is (f) or i R (f) is an even symmetric function. Then, the cross-correlation i{t) of is (t) and 
i R (t) is given as follows: 

i s (t) ® i R (0 = i s (0 * i R (0 = {e- 0 ^ 2 cos(? -t d )+w s (t) )*( e- fl ' 2 cos f + w ff (0 ) 
= jV al2 cos /I +w J? (A)) (e- a( '^ )2 cos(t-A-t d ) + w s (t-A))dA 

-00 

+co +oo (A- 12) 

= \((e~ al2 cos A) (e- fl (^') 2 costf-A-^))^ + jV^ cos 2 ) (w s (f - A) ) <M 

—00 -00 
+00 +00 

+ \(w R (A)) (e~ o(, "" A ~ , ' )2 cos(t-A-t d ))dA + \w R (A) w s (t-X))dA 



The first term in Equation A- 12 is the cross-correlation of noise free interference signal 
i NF (t) = e~ at cos(?) and its delayed version i NF (t-t d ) = e~ a( - l ~ td) cos(t-t d ), which is calculated as 
follows: 

i N F(t)®i NF (t-t d ) = i NF (t)*i NF (t-t d ) = {e- at2 cos^(e- fl( ^ )2 cos{t-t d )) (A-l 3) 
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Well-known Fourier transform relationships useful for the calculation of Equation A- 13 are: 

A(t)*B(t) <=> A(w)B(w) (A- 14) 

1 

2;r 



A(t)B(t)o—A(w)*B(w) (A- 15) 



7^4— ]»^(cw) (A-16) 
|«| \aJ 

where: 

^(0 o. ^(w) and B(t) «• 5(w) (A- 17) 

and w is the frequency domain variable. 

If we only consider the envelope of i NF {t) = e~ at cost , i.e. i NF env (t) = e~ at and use the Fourier 

transform relation: 

j w 2 

e - a ' 2 oJ^ (A-18) 
then Fourier transform of Equation A- 13 becomes: 

( \ I \ \ -— | -— I I (^w) 2 

^v(0*W, e „v(^^) = k" 2 H e_a( ^ )2 hd-e'\\-e~^)e-^ = (J- JV 4 « 

(A-19) 

| — I" (V2w) 2 

Then i NF env (t)*i NF env (t -t d )'vs the inverse Fourier transform of ( — — e 4a )e~ Jwt * , which is 

V a V a 

given as: 

Before cross-correlation, from Equation A- 10 and Equation A-ll, peak value of is (t) is 1 and is(t) 
reduces down to e 1 of its peak value at: 

t = ± A= (for Equation A- 1 0 ) and t = t d ± A= (for Equation A-ll) ( A-2 1 ) 



t = t d ±^= (A-22) 



After cross-correlation, from Equation A-20, e l of i s <8> i s (t) happens at: 

V2 

And the coherence length of the cross-correlation(effective coherence length) has increased by a 
factor of V2 . In other words, L c eff =42L C 

The second part of appendix A shows that 14 dB SNR improvement is obtained at zero order fringe 
peak after cross-correlation. Peak value (zero-order fringe peak) of cross-correlation i{t) happens when 
t = td and its value is* (t=td) is calculated by substituting t = td in Equation A- 12 as follows: 

+00 +co 

i s *i R (t d )= \((e~ al2 cos A) (e' a( - Af cos(-A))dA + \{e~ al1 cos A ) (w s (t d - A))dA 

(A-23) 

+00 +00 
+ \{w R {A)) (e- a( - A)1 cos{-A))dA + \w R (A) w s (t d -A))dA 
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The first term of Equation A-23 is the signal component which is the energy of noise free 
interference signal i NF (t) = e~ at cos(f) and the other three terms of Equation A-23 are noise terms. 

Then the variance of three noise terms are given as: 

+0O +00 

(2j(e^ 2 cos2) 2 dX)xa 2 N + (\w R (X) 2 dX)xa 2 N (A-24) 

—GO —00 

Then SNR at the peak of cross-correlation i (t) is given as the ratio between magnitude of peak 
signal to the standard deviation of noise as shown in Equation A-25. 

+CO 

\{e- aXl cos A) 2 d A 

SNR = , (A-25) 

I +GO +00 

2 jV fll2 cosX ) 2 dX) a 2 N +( \w R (X) 2 dX)a 2 N 

I -00 -00 

By assuming that the variance of the noise is much smaller than the power of signal: 

+00 +00 

\{e' al2 cos 2) 2 dX » \w R (X) 2 dX (A-26) 

-00 -00 

Equation A-25 can be approximated to: 



+00 



\{e al2 cosX) 2 dX \{e- al1 cosX) 2 dX 

SNR * , =g = = ( A-27) 

\{2\{e- al cos A ) dX)a\ 

\ -00 

Again, SNR of is (t) (or z«(0) before the cross-correlation is: 

snr Jnfmax_ (A _ 28) 

<J N 

so cross-correlation has improved the SNR by less than a factor of: 



+co 



I (e~ aA2 cos X ) 2 dX | e~ 2aXl cos 2 X dX 
SNR, lll/)mvc = ^ = x T ^= (A-29) 

"v 2 (Tjy l NF,MAX V2 l NF,MAX 

+00 

It can be shown that energy of the (noise-free) interference signal ^e 2aA cos 2 X dX in Equation 

-00 

A-29 is given by — — \=J— = — — J=J 7T 3 N 2 f as follows: 
2V2 \a 24l^ f 

00 - 00 

i NF* i N F (°)= \e~ 2axl cos 2 xdx =- \e- 2ax2 (\ + cos2x)dx (A-30) 

-co -00 

Useful relations are given by: 

1 [tt 



IV W coslxdx = —\-e 2a (A-32) 
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Then Equation A-30 becomes: 



—00 



+00 



cos 2 x dx = 



1 1 [tt 1 J_ l^r - 



2a 



(A-33) 



But, the second term in Equation A-29 is negligible because e 2a has a negligible value and 
Equation A-33 becomes: 



which is 14 dB when Nf = 26 
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